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ABSTRACT 


A  baseline  of  about  1200  km  between  Seattle,  Washington, 
and  Monterey,  California,  was  measured  repeatedly  over  a 
six-month  period  using  five-channel,  single-frequency  Global 
Positioning  System  (GPS)  receivers  with  carrier  phase 
differencing  techniques  and  broadcast  ephemeris.  The 
averaged  GPS  baseline  length  compared  favorably  with  the 
length  determined  from  control  points  established  by  Very 
Long  Baseline  Interferometry  (VLBI) ,  the  agreement  being  on 
the  order  of  0.01  ppm  (1  cm  in  1200  km)  which  is  about  the 
precision  expected  of  the  VLBI  technique  itself.  The 
quality  of  the  agreement  is  startling,  considering  the 
relatively  poorer  precision  (about  1  ppm)  expected  for  the 
GPS  receivers  and  techniques  employed.  To  achieve  this 
agreement,  GPS  observations  varying  more  than  1  ppm  from  the 
computed  mean  length  were  discarded,  and  a  scale  factor  of 
-0.2  ppm  for  the  transformation  from  GPS  to  VLBI  reference 
frames  was  applied,  which  had  been  estimated  from  other 
studies.  The  results  suggest  that  accuracies  of  better  than 
a  decimeter  are  achievable  over  lines  of  1000  km  using 
single-frequency  GPS  equipment. 
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I.   INTRODUCTION 

The  objective  of  this  study  is  to  determine  the  length 
of  a  baseline  between  Seattle,  Washington,  and  Monterey, 
California  using  two  independent  methods.  In  one  method, 
independent  Global  Positioning  System  (GPS)  surveys  were 
used  to  determine  Cartesian  coordinates  for  GPS  antennas  in 
Seattle  and  Monterey.  These  surveys  utilized  the  National 
Geodetic  Survey  (NGS)  Very  Long  Baseline  Interferometry 
(VLBI)  stations  located  near  the  Seattle  and  Monterey 
antenna  sites.  The  second  method  used  GPS  only  to  determine 
the  average  position  coordinates  and  slope  distance  between 
the  two  antennas  over  a  six-month  period. 

GPS  presently  consists  of  a  seven-satellite  constella- 
tion in  two  orbital  planes  separated  in  longitude  by  120° 
and  inclined  at  63°  to  the  eguatorial  plane.  On  average, 
three  to  six  satellites  are  in  view  from  two  to  six  hours 
each  day,  depending  on  geographical  location.  All  satel- 
lites are  at  an  altitude  of  20,183  km,  resulting  in  12-h 
sideral  orbits.  A  satellite  appears  in  view  approximately 
four  minutes  earlier  each  day  due  to  the  earth's  rotation  of 
almost  361°  every  24  hours. 

The  GPS  constellation  will  be  increased  to  18  satellites 
in  the  1990' s.  The  full  18-satellite  constellation  will 
make  it  possible  to  view  four  or  more  satellites  worldwide 


with  a  minimum  of  5°  elevation  above  the  horizon  (Stein, 
1986) . 

GPS  satellites  contain  cesium  and/or  rubidium  atomic 
frequency  standards  accurate  to  a  few  parts  in  1013  per  day 
(Stein,  1986) .  GPS  receivers  have  an  internal  quartz 
oscillator  clock.  The  receiver  beats  the  signal  from  its 
clock  with  the  signal  from  a  satellite  to  produce  a  beat 
frequency.  This  clock  comparison  is  called  the  carrier  beat 
phase  measurement. 

The  GPS  satellites  transmit  two  signals  simultaneously, 
LI  (1575.42  MHz)  and  L2  (1227.6  MHz).  LI  is  modulated  with 
a  precision  (P)  code  and  a  coarse  acquisition  (C/A)  code, 
whereas  the  L2  signal  is  modulated  with  either  the  P-  or 
C/A-code  (Stein,  1986)  .  The  P-  and  C/A-codes  are  used  to 
identify  the  GPS  satellites  and  for  transit  time  ranging 
between  the  receiver  and  the  satellites  (Remondi,  1984).  P- 
code  modulation  is  ten  times  faster  than  C/A-code  and 
provides  ten  times  greater  precision  (Smith,  1987) . 
Presently,  P-code  is  available  to  all  users  but  will  be 
limited  to  Department  of  Defense  once  the  network  is  fully 
operational . 

Mader  and  Abell  (1985)  found  baseline  lengths  to  be 
repeatable  to  within  an  average  of  0.24  ppm  for  distances  of 
300  to  1600  km  for  two-day  periods,  and  repeatabilities  of 
about  2  ppm  have  been  measured  over  300-  to  500-km  baselines 
over  five-day  periods  (Lachapelle  and  Cannon,  1986) . 


II.   METHODS 

A.   DIFFERENCING  THE  CARRIER  PHASE 

Correlations  may  be  made  among  signals  received  at 
stations  simultaneously  tracking  the  same  satellites. 
Errors  in  such  signals  may  also  be  correlated.  Because  of 
these  correlations  the  accuracy  of  relative  positions  may  be 
improved  by  taking  the  differences  between  measurements  to 
remove  or  greatly  reduce  errors. 

1.   Single  Difference 

Single  differencing  (Figure  1)  eliminates  receiver 
clock  errors.  A  single  difference  is  formed  by  differencing 
carrier  phases  at  the  same  epoch.  A  single  difference  can 
be  expressed  as  (Ashkenazi  et  al.,  1985): 

nABj(t)    =    ctABj    +    (fB   "    fA)(T  "    To)    +    f^CRfijCT) 

c 
AB  AB 

-    RAj(i)]    +    Nion    +    Ntrop  (1) 


where 


aABj    =    [NB(To)     "    NA(To)]     "    [NB,ro)     "    NA,ro)]  (2) 


Figure  1.   Single  Differencing  with  One  Satellite 
(  Wells  et  al. ,  1986  ) 


and 
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trop 


:   initial  (at  time  t0)  clock  bias; 

:   difference  at  local  epoch  t  of  the  GPS 

carrier  phase  readings  at  stations  A  and  B; 

:   freguency  offset  between  the  receivers  at 
stations  A  and  B; 

satellite  oscillator  freguency; 

initial  lock-on  time,  observation  time; 

velocity  of  microwave  propagation  in  vacuo; 

the  range  between  stations  A,B  to  satel]  I  Me 

J  1 

ionospheric  delay; 
tropospheric  delay. 


2 .   Double  Difference 

To  eliminate  both  local  and  satellite  clock  errors, 
double  differences  (Figure  2)  are  formed  by  differencing  :he 
single  differences  between  two  satellites  at  the  same  epoch. 
The  double  difference  can  be  written: 


Figure  2.   Double  Differencing  with  Two  Satellites 
( Wells  et  al. ,  1986  ) 


NABjk(T)     =    aABjk    +    NABjk    +    NABjk    +    f^R^,) 


"    RAkd)     ~    RBj(T>     +    RAj(^)l 


(3 


where 


NABjk('0     =    NABk(T)     "    NABj(T) 


(4 


and 


a 


ABjk(0  -  aABk(T)  "  «ABj(T) 


(5! 


3 .   Triple  Difference 

Triple  differencing  (Figure  3)  eliminates  all  errors 
found  in  double  differencing  as  well  as  cycle  slips.  The 
triple  difference  can  be  written: 


NABjk(Tl'T2)  =  NABjk(T2)  "  NABjk(Tl) 


(e 


Figure  3.   Triple  Differencing  with  Four  Satellite 
(  Wells  et  al. ,  1986  ) 


B.   ERRORS  IN  SINGLE-FREQUENCY  DATA 

The  sum  of  all  biases  adding  to  or  subtracting  frcm  a 
satellite's  projected  range  is  called  the  range  bias. 
Common  biases  found  in  single-frequency  observations  are 


related  to  satellite  clock,  receiver  clock,  orbital,  iono- 
spheric and  tropospheric  delays,  and  the  carrier  beat  phase. 

Satellite  clock  biases  may  lead  to  10-m  range  errors  if 
broadcasted  corrections  are  used,  while  receiver  clock 
biases  may  lead  to  10-  to  100-m  errors,  depending  on  the 
type  of  receiver  oscillator  (Wells  et  al.,  1986).  As 
explained  above,  these  errors  can  be  removed  by  differencing 
the  data. 

Orbital  biases  result  from  the  departure  of  a  satellite 
from  its  broadcasted  ephemeris  or  predicted  orbit.  Orbital 
biases  propagate  into  a  computed  baseline  when  the  GPS  orbit 
coordinates  are  fixed  in  processing  (Hothem  and  Williams, 
1985) ,  and  can  lead  to  errors  of  up  to  80  m  for  broadcasted 
ephemerides  (Wells  et  al.,  1986). 

The  GPS  signal  is  affected  by  nonlinear  dispersion  in 
the  atmosphere  which  can  lead  to  range  biases  of  more  than 
150  m  at  a  sunspot  maximum  to  less  than  5  m  at  a  sunspot 
minimum  (Wells  et  al.,  1986).  Group  velocities  at  radio 
frequencies  are  retarded  by  the  ionosphere  by  an  amount 
proportional  to  the  total  electron  content  (TEC)  along  the 
signal's  path  and  inversely  proportional  to  the  square  of 
the  frequency  (Goad,  1985)  .  The  TEC  is  a  function  of  the 
time  of  day  and  year,  latitude,  longitude,  and  sunspot 
activity  (Stoichar,  1985)  .  Ionospheric  errors  increase 
towards  the  equator  where  sunlight  is  more  intense.  Dual- 
frequency  receivers  can  be  used  to  measure  the  ionospheric 


delay  by  comparing  delays  at  the  two  frequencies  (Kaniuth, 
1986) .  For  single-frequency  receivers  ionospheric  delay 
must  be  estimated  on  the  basis  of  TEC  predictions  (Smith, 
1987)  . 

Tropospheric  biases  are  proportional  to  the  refractivity 
found  in  the  non-ionized  atmosphere  along  the  satellite- 
receiver  path.  Radio  waves  in  the  stratosphere  and  tropo- 
sphere are  not  dispersive  up  to  30  GHz  (Kaniuth,  1986)  . 
Refractivity  can  be  written  as  N  =  (n-1)  x  106  where  n  is 
the  refractive  index.  Tropospheric  biases  vary  from 
approximately  2.3  m  at  the  zenith  to  2  0  m  at  10°  above  the 
horizon. 

Carrier  beat  phase  biases  can  result  in  gross  errors.  A 
1-ys  miss-synchronization  between  the  satellite  and  receiver 
clocks  creates  a  300-m  range  bias  (Wells  et  al.,  1986). 

Position  Dilution  of  Precision  (PDOP)  is  a  component  of 
Geometric  Dilution  of  Precision  (GDOP) ,  which  is  a  measure 
of  how  satellite  geometry  degrades  accuracy  (Jorgenson, 
1984).  The  Trimble  4000SX  receiver  used  for  our  observa- 
tions records  PDOP  every  five  minutes.  PDOP  is  related  to 
GDOP  by: 


GDOP2  =  PDOP2  +  TDOP2  (7) 


where  TDOP  is  Time  Dilution  of  Precision  (Jorgenson,  1984) . 
TDOP  is  the  error  in  the  user  clock  bias  multiplied  by  the 


propagation  speed.   PDOP  peaks  occur  when  the  satellites  lie 
in  a  common  plane. 

GPS  observations  should  be  made  close  to  2400  local  time 
once  24-h  satellite  visibility  is  available.  This,  in 
conjunction  with  low  PDOPs,  should  minimize  ionospheric  and 
tropospheric  errors. 


III.   SITE  SURVEYS  AND  INSTRUMENTS 

A.   SEATTLE,  WASHINGTON 

The  first  factors  to  determine  in  establishing  the 
baseline  were  the  positions  of  the  antennas  in  Seattle  and 
Monterey.  This  was  done  by  using  Trimble  4000SX  GPS 
receivers.  One  receiver  was  installed  and  tested  on  25 
September  1987  at  NOAA's  Sand  Point  Facility  in  Seattle, 
Washington  (Figures  5  and  6)  .  A  second  Trimble  4000GPS 
receiver  was  used  to  locate  the  Sand  Point  antenna  with 
respect  to  a  nearby  VLBI  control  point  (referred  to  as 
AVIATION  2).  This  VLBI  antenna  position  and  another  near 
the  Monterey  GPS  antenna  site  were  used  to  check  the  GPS- 
determined  baseline  between  Seattle  and  Monterey  (Figure  4). 
Satellites  6,  9,  11,  12,  and  13  were  used  during  our 
observations. 

The  GPS  antenna  position  at  Seattle  was  determined  using 
one  90-min  GPS  observation  session  with  a  double  difference 
carrier  phase  solution  (Table  1) .  This  solution  determines 
the  change  in  coordinates  between  the  two  antennas.  Let  A 
be  the  difference  between  the  corresponding  coordinates  for 
the  antennas  at  stations  Aviation  2  and  Seattle  and  -,  its 
standard  deviation.  The  standard  deviation  of  Seattle,  eg, 
is : 
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Figure  4.   Station  Location  Map  Showing  Monterey, 
California,  and  Seattle,  Washington 
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Figure  5.   The  Environs  of  Seattle,  Washington,  Showing 
the  Location  of  Sand  Point  and  the  Seattle- 
Tacoma  Airport 
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Figure  6.   Location  of  Seattle  GPS  Antenna  Relative 
to  the  VLBI  Control  Point 


TABLE  1 

RESULTS  OF  ANTENNA  LOCATION  SURVEY  IN  SEATTLE 

(Bouchard,  1988) 


Aviation  2 
Coordinate 

X  -2295347.760 
Y  -3638029.429 
Z    4693408.964 


°A 
0.017 
0.028 
0.032 


Seattle 
Coordinate 

-2295756.121 

-3637699.228 

4693482.777 


US 
0.017 
0.028 
0.032 


-408.361 

330.200 

73.813 


0.002 
0.  002 
0.003 


Receiver  clock,   satellite  clock,   and  orbital  errors  are 
likely  responsible  for  the  standard  deviations  in  Table  1. 
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B.   MONTEREY,  CALIFORNIA 

The  GPS  station  position  at  Monterey  (Figure  7)   was 
found  by  differential  GPS   using   FT   ORD  NCMN  1981,  a  VLBI 
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Figure  7.   Location  of  Monterey  GPS  Antenna  Relative 
to  the  VLBI  Control  Point 
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control  point,  as  the  primary  station.  Data  collected 
during  two  observing  sessions  were  averaged  to  give  antenna 
positions.  The  first  observing  session  was  112  min  in 
length  while  the  second  was  126  min.  The  Cartesian 
coordinates  of  each  GPS  antenna  used  in  the  observing 
sessions  were  included  in  the  post-processed  output.  These 
values  are  listed  in  Table  2  under  the  columns  FT  ORD  and 
Monterey.  a  and  oA  were  computed  as  described  in  Section  A 
above. 

TABLE  2 
RESULTS  OF  ANTENNA  LOCATION  SURVEYS  IN  MONTEREY 


SURVEY  DATE:   09/16/87 

FT  ORD       a0  Monterey  aM  A  aA 

X  -2697026.493   0.007  -2707340.093  0.032  -10313.601  0.032 

Y  -4354393.309   0.010  -4353475.617  0.049  917.693  0.048 

Z   3788077.778   0.009  3781740.387  0.042  -6337.391  0.041 


SURVEY  DATE:   09/18/87 


FT  ORD 

aQ 

X 

-2697026.493 

0.007 

Y 

-4354393.309 

0.010 

Z 

3788077.778 

0.009 

Monterey  aM  A  o^ 

-2707340.095  0.037  -10313.602  0.037 

-4353475.620  0.050  917.689  0.049 

3781740.389  0.045  -6337.389  0.044 


AVERAGE  MONTEREY  ANTENNA  POSITION: 

Monterey  o^ 

X  -2707340.094  0.034 

Y  -4353475.618  0.050 

Z   3781740.388  0.044 
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IV.   DATA  ACQUISITION 

A.   GPS 

GPS  positioning  data  acquired  throughout  the  experiment 
were  collected  using  at  least  four  of  the  five  satellites 
available  (6,  9,  11,  12,  and  13).  This  set  was  chosen 
because  the  satellites  were  well  positioned  for  viewing  from 
both  Seattle  and  Monterey  (Figure  8)  .  Five  satellites 
allowed  observation  periods  of  at  least  100  min.  For  an 
additional  80  min  four  satellites  were  still  visible. 

HOIUII 


Figure  8.  Satellite  Tracks  Relative  to  Monterey.  The 
Zenith  for  Monterey  is  at  the  Center,  While 
the  Horizon  is  on  the  Outer  Perimeter 
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GPS  data  were  collected  from  29  September  1987  until  31 

March  1988  with  Trimble  4000SX  single-frequency  (C/A-code) , 

five-channel  receivers.   The  receivers  log  all  carrier  phase 

data  via  Grid  laptop  microcomputers  onto  3.5-in  floppy  disks 

(Trimble  Navigation,  1987) . 

Observations  were  made  Tuesdays  through  Saturdays  except 
on  days  after  Federal  holidays.  Reasons  for  not  processing 
certain  days  included  bad  satellite  health,  power  failures, 
or  the  unavailability  of  meteorological  data. 

The  4  000SX  receiver  uses  the  C/A  code  to  decode  the  GPS 
navigation  message  so  it  can  automatically  track  satellites. 
It  also  acts  in  a  time  transfer  mode  to  determine  offset  and 
drift  of  its  own  clock,  thus  providing  accurate  time  tags 
for  observations  without  external  atomic  clocks  or  synchron- 
ization with  the  receiver  at  the  other  end  of  the  baseline. 
The  receivers  were  left  on  continuously  throughout  the  six- 
month  period  to  allow  unattended  data  acquisition.  They 
were  controlled  using  Trimble's  Version  D  "Datalogger" 
acquisition  software.  The  station  reference  position, 
antenna  height,  and  any  additional  information  were  entered 
into  the  receiver  via  its  keypad. 

Each  day  the  desired  satellites  and  their  minimum 
acceptable  elevations  were  set  from  the  computer.  A  minimum 
elevation  of  15  degrees  above  the  horizon  was  used  to 
control  when  the  receivers  would  start  collecting  satellite 
data.    This  15-degree  elevation  is  based  on  when  minimum 
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error  biases  due  to  both  the  ionized  and  non-ionized  parts 
of  the  atmosphere  start  to  occur. 

Receiver  clock  parameters  were  logged  to  the  floppy  disk 
every  15  s  and  antenna  positions  every  5  min.  The  GPS 
navigation  message,  ionospheric  and  Universal  Coordinated 
Time  (ION/UTC)  data,  and  an  optional  user  message  were 
transferred  to  separate  files  at  the  beginning  of  each 
session.  The  ION/UTC  data  contain  ionospheric  parameters 
for  dual  frequency  receivers  and  UTC  time  parameters  to 
initialize  the  receiver's  clock. 

B.   METEOROLOGICAL 

Meteorological  parameters  for  tropospheric  corrections 
were  obtained  daily  from  the  NPS  Department  of  Meteorology. 
Values  for  Seattle  were  recorded  by  the  Weather  Service 
Office  of  the  Seattle-Tacoma  International  Airport  (SEA-TAC) 
approximately  27  km  from  the  Seattle  site  (Figure  5)  . 
Values  used  for  Monterey  were  measured  hourly  by  the  NPS 
Department  of  Meteorology  approximately  300  m  south  of  the 
antenna  site. 
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V.   DATA  ANALYSIS 

A.   COORDINATE  SYSTEM 

All  GPS  data  were  processed  using  the  World  Geodetic 
System  1984  (WGS  84)  coordinate  system,  based  upon  broad- 
cast ephemeris  of  the  GPS  satellites.  However,  a  small  but 
systematic  difference  between  GPS-  and  VLBI-derived  WGS  84 
coordinates  turned  out  to  be  significant.  The  broadcast  GPS 
ephemeris,  which  predicts  positions  of  the  satellites  in  WGS 
84  coordinates,  affected  the  baseline  distances  determined 
by  the  carrier  phase  differential  technique.  On  the  other 
hand,  VLBI  control  points,  originally  in  the  VLBI  reference 
frame,  were  used  to  derive  WGS  84  coordinates  for  the  GPS 
antenna  positions  in  Seattle  and  Monterey,  and  hence 
affected  the  baseline  distance  estimated  by  differencing 
these  positions. 

The  possibility  of  systematic  differences  between  the 
GPS-  and  VLBI-derived  WGS  84  positions  has  been  investigated 
by  Abusali  et  al.  (1989).  In  their  study,  a  set  of 
positions  expressed  in  the  two  coordinate  systems  (VLBI  and 
GPS/WGS-84)  were  compared  to  empirically  determine  a  best- 
fitting  transformation  between  them.  A  conventional  seven- 
parameter  transformation  model  was  used:  three  translations 
and  three  rotations  of  the  XYZ  axes  plus  a  scale  factor  for 
any  proportional  differences  in  length.   The  specific  set  of 
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positions  used  happened  to  be  those  of  the  GPS  satellites 
determined  by  simultaneous  tracking  from  two  sets  of  ground 
stations,  one  with  coordinates  given  in  the  GPS  system  and 
the  other  given  in  VLBI-derived  coordinates. 

Their  results  suggested  that  a  scale  factor  on  the  order 
of  -0.1  to  -0.2  ppm  may  exist  (in  the  sense  of  transforming 
from  distances  in  the  GPS  system  to  those  in  the  VLBI 
system) .  However,  these  results  were  not  applicable 
directly  to  the  question  here;  that  is,  what  if  any  syste- 
matic scale  factor  exists  between  the  broadcast  GPS 
ephemeris  and  VLBI  systems?  The  broadcast  ephemeris 
includes  an  orbit  prediction  model  which  was  not  considered 
in  the  Abusali  study.  In  a  private  communication,  Abusali 
(1989a)  indicated  a  scale  factor  of  about  -0.2  ppm  (GPS  to 
VLBI^  should  be  expected  if  using  broadcast  ephemeris. 

This  estimate  of  the  scale  factor  was  applied  in  the 
comparison  between  GPS-  and  VLBI-determined  slope  distances 
as  discussed  in  Section  VI,  Results.  Otherwise,  it  was  not 
applied  to  the  data  analysis  described  here. 

B.   GPS-BASELINE  DETERMINATION 

1 .   Six-month  Averaged  Baseline  Length 

Observations  collected  during  each  GPS  session  were 
divided  into  10-min  segments  using  INTER_H.BAS  (Appendix  A) 
to  observe  the  effect  of  varying  satellite  geometry  on  the 
baseline.  GETHDAT.BAS  (Appendix  B)  was  used  to  read  and 
store  all  slope  distances  processed  by  Trimble's  TRIM640.EXE 
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(Revision  AB)  .    Outlying  values  of  slope  distance  were 
removed   by   eliminating   those   exceeding   a  3c   value.   3  o 
represents  99%  of  those  accepted  slope  distances  surrounding 
the  mean.   The  30  value  was  computed  by: 

EVi2  1/2 
where: 


Vj_  =  v  -   Xi; 

y  =  mean  value  of  slope  distance; 

Xj_  =  observed  value  of  slope  distance  based  on 
10-min  averages; 

n  =  total  number  of  values. 


Once  the  value  of  o  is  known,  it  can  be  found  in  parts  per 
million  (ppm)  as: 

(c/L)106  (10) 

where: 

L  =  the  mean  slope  distance. 

The  rejection  criterion  above  was  repeated  using  those 
values  within  3o  if  the  value  of  o  was  more  than  1  ppm. 
This  value  was  chosen  because  it  was  well  within  the 
standards  for  first  order  surveys  of  baselines  of  this 
length.    (For  additional  comments  see  Bouchard  (1988  ,  p. 
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50).)  Four  iterations  were  required  to  achieve  an  accuracy 
better  than  1  ppm.  The  standard  deviation  around  the  mean 
was  then  computed  by: 

°m  =  o/(n)V2  (11) 

Using  o m  in  place  of  o  in  Equation  (10),  ® m  may  be  expressed 
in  ppm.  The  final  computed  slope  distance  is  1230045.532 
±  0. 025  m  (0. 02  ppm)  . 

This  value  is  well  within  the  allowable  limit  for 
first  order  surveys  (1:100,000)  of  the  National  Geodetic 
Reference  System  (Terrestrial) . 

2 .   Ninety-Four-Day  Baseline  Length 

GPS  data  from  94  days  were  processed  using  triple 
differencing.  These  data  were  not  divided  as  in  Section  I 
but  rather  processed  at  one  time  to  ensure  that  the  PDOP 
passed  through  infinity.  It  is  with  this  rapidly  changing 
PDOP  that  best  results  are  acquired  (Trimble  Navigation, 
1987;  Bouchard,  1988).  The  average  X,  Y,  and  Z  values  are 
noted  in  Table  3.   The  slope  distance  was  computed  using: 

D  =  (AX2  +  AY2  +  AZ2)1/2  (12) 
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TABLE  3 
CARTESIAN  COORDINATES  OF  GPS  DETERMINED  ANTENNA  POSITIONS 


Seattle 
X      -2295759.503  ±  1.277 
Y      -3637704.115  ±  0.971 
Z       4693484.159  ±  0.729 


Monterey 
-2707343.469  ±  1.341 
-4353480.402  ±  0.945 
3781742.430  ±  0.731 


The  standard  error  in  the  slope  distance  was  then  computed 
as: 


D  =  U(AX  °AX)/D]2  +  t(AY  oAY)/D] 

[(AZa.z)/D]2}V2 


(13) 


where  D  is  the  slope  distance.    °AX'  ctaY'   anc^  °LZ     were 
computed  using: 


=  ±  (- ±2    +   o22) 1/2 


(14) 


where  c1  and  c2  are  the  standard  errors  of  Seattle  and 
Monterey. 

The  final  computed  slope  distance  was  1230044.750 

±  1.261  m.    It  is  suspected  that  the  high  PDOP  values 

included  in  the  data  processing  contributed  to  a  higher 
standard  deviation  for  D  when  compared  to  the  slope  distance 
in  Section  I. 
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3 .   VLBI-Determined  Antenna  Position 

The  VLBI  positions  were  obtained  in  NAD83  Cartesian 
coordinates  from  the  Gravity,  Astronomy  and  Space  Geodesy 
Branch  of  the  National  Geodetic  Survey  (NGS) .  The  Defense 
Mapping  Agency  Hydrographic/Topographic  Office  validated  the 
direct  transformation  of  the  VLBI  Cartesian  coordinates  to 
WGS84  Cartesian  coordinates  (Kumar,  1988)  . 

The  Monterey-Seattle  VLBI  antenna  positions  given  in 
Tables  1  and  2  were  used  to  check  the  baselines  computed  in 
Sections  I  and  II.  The  Cartesian  coordinates  of  the  VLBI- 
determined  antenna  positions  were  computed  using  double 
difference  processing.  This  processing  technigue  was  chosen 
because  the  root-mean-sguare  (RMS)  of  fit  was  less  than  0.05 
cycles,  which  is  common  for  baselines  less  than  3  0  km 
(Trimble  Navigation,  1987) .  The  method  used  to  compute  the 
Monterey-Seattle  baseline  and  its  standard  deviation  is 
discussed  in  Section  II.   Computed  values  are  in  Table  4. 

TABLE  4 
VLBI/GPS  POSITIONED  ANTENNA  BASELINE  DISTANCES 

AX  -411583.973  ±  0.040 

AY  -715776.390  -     0.058 

LZ  -911742.338  =  0.054 

Slope  Distance:      1230045.280    0.054 
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This  slope  distance  yields  0.04  ppm.  Precisions  of 
0.01  ppm  have  been  achieved  with  VLBI  and  1  ppm  with  GPS 
over  1600-km  baselines  in  Alaska  and  Canada  (Mader  and 
Abell,  1985) . 
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VI.   RESULTS  AND  CONCLUSIONS 

The  six-month  averaged  GPS  slope  distance  differs 
significantly  from  that  determined  by  the  VLBI-derived 
antenna  positions,  the  GPS  distance  being  about  25  cm  too 
long. 

Slope  Distance         Estimated 

(meters)  Standard  Deviation 

VLBI  1230045.280  0.054 

GPS  1230045.532  0.025 

Difference  -0.252 

Given  the  estimated  precisions  on  each  distance,  this 
difference  suggests  that  some  systematic  effects  remain. 
(The  hypothesis  of  equal  distances  can  be  rejected  at  about 
the  98%  confidence  level  in  favor  of  the  GPS  value  being 
greater  by  using  a  Fischer-Barents  test  for  samples  with 
unequal  variances.) 

As  discussed  in  Section  IV. A,  studies  at  the  Center  for 
Space  Research,  University  of  Texas  (Abusali,  et  al . ,  1989), 
suggest  that  a  systematic  scale  error  exists  between  the 
GPS-  and  VLBI-derived  coordinate  systems.  Adopting  their 
best  estimate  of  -0.2  ppm  for  the  scale  error  between  GPS 
broadcast  ephemeris  and  VLBI  coordinates  (Abusali,  1989), 
brings  these  measurements  into  close  agreement: 
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GPS  length   =   1230045.532  m 

-0.2  ppm   =   -0.  246  m 

1230045.286  m  Transformed  GPS  distance 

which  differs  by  about  1  cm  from  the  VLBI-derived  distance 
(with  an  estimated  precision  of  about  6  cm,  one  standard 
deviation) . 

This  result  indicates  that  accuracies  of  better  than  10 
centimeters  can  be  achieved  with  single-frequency  GPS 
equipment  over  relatively  long  baselines,  on  the  order  of 
1000  km.  Careful  averaging  of  repeated  measurements  is 
essential  in  order  to  reduce  the  incompletely  modelled 
effects  of  ionospheric  and  tropospheric  paths  differences 
over  such  long  lines.  And  although  this  result  was  achieved 
using  broadcast  ephemeris,  precise  ephemeris  should  be  used 
whenever  possible  in  order  to  eliminate  systematic  errors 
due  to  orbit  predictions  methods  employed  by  GPS. 

This  result  also  tends  to  confirm  the  existence  of  a 
systematic  scale  error  between  the  GPS  and  VLBI  coordinate 
systems.  This  scale  error,  however,  is  significant  only  to 
those  users  mixing  VLBI-  and  GPS-derived  information.  Most 
will  work  entirely  within  the  WGS  84  coordinate  system  using 
only  GPS-derived  positions. 
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APPENDIX  A 
PROGRAM  INTER  H . BAS  LISTING 


1  REM  THIS  PROGRAM  CREATES  A  BATCH  FILE  TO  BE  USED  WITH 
TRIMBLES  POSTPROCESSING  SOFTWARE  "TRIMVEC."   IT  DOES  A 
LINEAR  INTERPOLATION  OF  METEOROLOGICAL  VALUES  FOR  A 
GIVEN  DAY  FOUND  IN  THE  "MET. DAT"  FILE.   THIS  PROGRAM 
COMPUTES  FOR  EVERY  10  MINUTES. 

10  INPUT  "BATCH  FILE  NAME?";A$ 

2  0  OPEN  "DEPOSIT.H"  FOR  APPEND  AS  #1 

3  0  OPEN  "MET. DAT"  FOR  INPUT  AS  #2 

4  0  INPUT  "MONTH?" ;M$ 
50  INPUT  "DAY?";D$ 

60  INPUT  "JULIAN  DAY?";JD 

62  DIM  A$(18) 

63  A$(2)="  h02.":A$(3)="  h03 . " : A$ (4 ) ="  h04 . " : A$ (5) =" 

h05.":A$(6)="  h06 . " : A$ (7) ="  h07 . " : A$ (8) =" 

h08.":A$(9)=" 

h09.":A$(10)="  hlO.":A$(ll)="  hll. " : A$ (12) =" 

hl2.":A$(13)=" 

hl3.":A$(14)="  hl4.":A$(15)="  hl5 . " : A$ (16) =" 

hl6.":A$(17)=" 

hl7.":A$(18)="  hl8." 
65  JD$=STR$(JD) 
67  W$="  " 

70  C$="command  /c  tbf  h.tem  " : F$="sa"+JD$+" 
ma"+JD$:S$=C$+F$ 
80  INPUT  "START  HOUR?";SH 
90  INPUT  "START  MINUTE?" ;SM 
100  INPUT  #2,ID1,IH1,P1,T1,R1,P2,T2,R2 
110  IF  (IDloJD)  THEN  100 
120  IF  (SHOIH1)  THEN  100 
130  Pl=Pl-2:T2=(T2-32) *5/9 
140  INPUT  #2,ID2/IH2,P3,T3,R3,P4,T4,R4 
150  P3=P3-2:T4=(T4-32)*5/9 
160  IF  ID10ID2  THEN  260 
165  REM  THIS  NEXT  SECTION  COMPUTES  FRACTIONAL  VALUES 

FOR  FUTURE  LINEAR  INTERPOLATION. 
170  FOR  1=1  TO  6 

180  A=A+IH1 : B=B+P1 : C=C+T1 : D=D+R1 : E=E+P2 : F=F+T2 : G=G+R2 
190  PRINT  #1,ID1,A,B,C, D,E,F,G 
200  PRINT  ID1,A,B,C,  D,E,F,G 
210  IF  I>1  GOTO  230 
220 
IHl=(IH2-IHl)/6:Pl=(P3-Pl)/6:Tl=(T3-Tl)/6:Rl=(R3-Rl)/6: 

P2=(P4-P2)/6:T2=(T4-T2)/6:R2=(R4-R2)/6 
2  30  NEXT  I 
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24  0  ID1=ID2 : IH1=IH2 : P1=P3 : T1=T3 : R1=R3 : P2=P4 : T2=T4 : R2=R4 

245  A=0:B=0:C=0:D=0:E=0:F=0:G=0 

250  GOTO  140 

260  CLOSE  #1 

270  CLOSE  #2 

280  OPEN  A$  FOR  APPEND  AS  #1 

290  OPEN  "DEPOSIT. H"  FOR  INPUT  AS  #2 

300  INPUT  #2,ID1,IH1,P1,T1,R1,P2,T2,R2 

310  IF  (IDloJD)  THEN  300 

320  IF  (SHOINT(IHI)  )  THEN  300 

3  25  REM  FRACTIONAL  MINUTE  IS  NOW  COMPUTED  AND  TESTED 

WITH  THE  FRACTIONAL  HOUR. 
330  FM=SM/60 

340  IF  (SH+FM-.0833)>IH1  OR  (SH+FM+ . 0833 ) <IH1  THEN  300 
344  IH=INT(IH1) : IM= (IH1-IH) *60 

346  IF  IM-INT(IM)>.5  THEN  IM=IM+ ( 1- (IM-INT (IM) ) ) 
348  IF  IM-INT(IM)<.5  THEN  IM=INT(IM) 

350  IM2=IM+10:IH2=IH 

351  REM  THIS  NEXT  SECTION  IS  A  ROUNDING  OFF  ROUTINE. 

352  IF  IM2>60  THEN  IH2=IH+1  AND  IM2=IM2-60 
354  IF  P1=INT(P1)  THEN  356 

3  55  Pl=INT(Pl)+( (INT( (Pl-INT(Pl) )*10) )/10) 

356  IF  T1=INT(T1)  THEN  358 

3  57  Tl=INT(Tl)+( (INT( (Tl-INT(Tl) )*10) )/10) 

358  R1=INT(R1) 

359  IF  P2=INT(P2)  THEN  361 

3  60  P2=INT(P2)+( (INT( (P2-INT(P2) ) *10) )/10) 

361  IF  T2=INT(T2)  THEN  363 

3  62  T2=INT(T2)+( (INT( (T2-INT(T2) ) *10) )/10) 

363  R2=INT(R2) 

370 

P1$=STR$(P1) :T1$=STR$(T1) : R1$=STR$ (Rl) : P2$=STR$ (P2 ) : 

T2$=STR$(T2) :R2$=STR$(R2) : P$="h01 . " : IM$=STR$ (IM) :I- 

H$=STR$(I- 

H) :IM2$=STR$(IM2) : IH2$=STR$ (IH2) 
380  PRINT 

#1,S$+P$+JD$+P1$+T1$+R1$+P2$+T2$+R2$+W$+M$+W$+D$+- 

IH$+IM$+IH2$+IM2$ 
390  PRINT 

S$+P$+JD$+P1$+T1$+R1$+P2$+T2$+R2$+W$+M$+W$+D$+ 

IH$+IM$+IH2$+IM2$ 
400  IM=IM2:IH=IH2 
410  FOR  N=2  TO  18 
420  INPUT  #2,ID1,IH1,P1,T1,R1,P2,T2,R2 

423  IH=INT(IH1) : IM= (IH1-IH) *60 

424  IF  IM-INT(IM)>.5  THEN  IM=IM+ ( 1- (IM-INT (IM) ) ) 

425  IF  IM-INT(IM) <.5  THEN  IM=INT(IM) 

426  IM2=IM+10:IH2=IH 

427  REM  THIS  NEXT  SECTION  IS  A  ROUNDING  OFF  ROUTINE. 

429  IF  IM2>60  THEN  IH2=IH+1  AND  IM2=IM2-60 

430  IF  P1=INT(P1)  THEN  432 

431  Pl=INT(Pl)+( (INT( (Pl-INT(Pl) ) *10) )/10) 
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432  IF  T1=INT(T1)  THEN  434 

433  Tl=INT(Tl)+( (INT( (Tl-INT(Tl) )*10) )/10) 

434  R1=INT(R1) 

435  IF  P2=INT(P2)  THEN  437 

4  36  P2=INT(P2)+( (INT( (P2-INT(P2) ) *10) )/10) 

437  IF  T2=INT(T2)  THEN  439 

4  38  T2=INT(T2)+( (INT( (T2-INT(T2) )*10) )/10) 

439  R2=INT(R2) 

447 

P1$=STR$(P1) :T1$=STR$(T1) :R1$=STR$ (Rl) : P2$=STR$ (P2) : 

T2$=STR$(T2) :R2$=STR$(R2) : IM$=STR$ (IM) : IH$=STR$ (IH- 

) :IM2$=ST-R$(IM2) : IH2$=STR$ (IH2) 
450  PRINT 

#1,S$+A$(N)+JD$+P1$+T1$+R1$+P2$+T2$+R2$+W$+M$+ 

W$+D$+IH$+IM$+IH2$+IM2$ 
460  PRINT  S$+A$(N)+JD$+P1$+T1$+R1$+P2$+T2$+R2$+W$+M$+ 

W$+D$+IH$+IM$+IH2$+IM2$ 
470  IM=IM2:IH=IH2 
480  NEXT  N 
490  END 
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APPENDIX  B 
PROGRAM  GETHDAT.BAS  LISTING 


1  REM  THIS  PROGRAM  WITHDRAWS  THE  SLOPE  DISTANCE,  X,  Y, 

AND  Z,  OF  TRIMVEC'S  ASCII  PRINTOUT. 

4  DATA  272,273,274,275,276,279,280,281,282,287 

5  FOR  K  =  1  TO  46 
7  A  =  20 

10  READ  JD$ 

2  0  OPEN 

2  5  OPEN 

3  0  OPEN 

3  5  OPEN 

4  0  OPEN 
4  5  OPEN 
50  OPEN 
55  OPEN 
60  OPEN 
65  OPEN 
7  0  OPEN 

7  5  OPEN 

8  0  OPEN 

8  5  OPEN 

9  0  OPEN 
9  5  OPEN 


"HOI." 

,,H02.,, 

"H03. ' 

"H04.  ' 

"H05.' 

"H06.1 

"H07.1 

"H08.1 

"H09.' 

"H10.' 

"Hll." 

"H12.1 

"H13.' 

"H14. ' 

"H15.' 

"H16.' 

•+JD$  FOR    INPUT    AS 

'+JD$  FOR    INPUT   AS 

•+JD$  FOR    INPUT   AS    #1 

■+JD$  FOR    INPUT   AS    #1 

'+JD$  FOR    INPUT   AS 

•+JD$  FOR    INPUT   AS 


#1 
#1 


#1 
#1 


!+JD$  FOR  INPUT  AS  #1 
•+JD$  FOR  INPUT  AS  #1 
'+JD$  FOR  INPUT  AS  #1 


GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 
GOTO  110 


FOR  APPEND  AS  #2 


'+JD$  FOR  INPUT  AS  #1 
•+JD$  FOR  INPUT  AS  #1 
•+JD$  FOR  INPUT  AS  #1 
■+JD$  FOR  INPUT  AS  #1 
'+JD$  FOR  INPUT  AS  #1 
•+JD$  FOR  INPUT  AS  #1 
'+JD$  FOR  INPUT  AS  #1 

100  OPEN  "H17."+JD$  FOR  INPUT  AS  #l:GOTO  110 

105  OPEN  "H18.U+JD$  FOR  INPUT  AS  #l:GOTO  110 

110  A=A+5 

115  OPEN  "GETHDAT.PRN" 

12  0  FOR  CA  =  1  TO  14 

122  INPUT  #1,  A$ 

12  4  NEXT  CA 

12  6  INPUT  #1,  A$ 

128  W$=MID$(A$,12,10) 

129  IF  A  >=  95  GOTO  144 

13  0  FOR  CA  =  1  TO  7  8 
132  INPUT  #1,  A$ 

14  0  NEXT  CA 
142  GOTO  150 
14  4  FOR  CA  =  1  TO  7  4 
146  INPUT  #1,  A$ 
14  8  NEXT  CA 
150  FOR  CJ  =  1  TO  4 
160  INPUT  #1,  A$ 
170  IF  CJ  =  2  GOTO  200 
180  IF  CJ  =  3  GOTO  210 
185  IF  CJ  =  4  GOTO  215 
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190 
200 
210 
215 
220 
240 
260 
270 
272 
280 
300 
310 
320 
330 
335 
340 
345 
350 
355 
360 
365 
370 
375 
380 
385 
390 
395 
400 
405 
410 
415 


AA$  =  MID$(A$, 11,12) 
B$  =  MID$(A$,11,12) 
C$  =  MID$(A$, 11,12) 
H$  =  MID$(A$,18,8) 
NEXT  CJ 
INPUT  #1,  A$ 
INPUT  #1,  A$ 
G$  =  MID$(A$,24,12) 
PRINT  W$,  G$,  C$,  H$ 


PRINT  #2, 
CLOSE  #2 
CLOSE  #1 


W$,  G$,  C$,  H$ 


IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

NEXT 

END 


25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 


GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 


25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 


=  105  GOTO  105 
K 
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